perm filename SPECT.F4[1,MUS] blob sn#007307 filedate 1972-07-16 generic text, type T, neo UTF8
00100		DIMENSION F(20)
00200		DO 1 I=1,20
00300	1	F(I)=0
00400		ACCEPT 100,CF,XMF
00500	100	FORMAT (2F)
00600		F(1)=CF
00700		J=2
00800		DO 2 I=1,9
00850		X=I
00900		F(J)=CF+(X*XMF)
01000		J=J+1
01100		F(J)=CF-(X*XMF)
01200		J=J+1
01225	2	CONTINUE
01250		DO 3 K=1,18
01300		DO 3 I=1,18
01400		X=ABS(F(I))
01500		X2=ABS(F(I+1))
01600		IF(X.LT.X2)GO TO 3
01700		X=F(I)
01800		F(I)=F(I+1)
01900		F(I+1)=X
02000	3	CONTINUE
02050		TYPE 102,CF,XMF
02075	102	FORMAT(' CF='F8.3,4X'MF='F8.3/)
02100		TYPE 101,F
02200	101	FORMAT(F10.3)
02300		END